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Abstract 

We present dynamical properties of spin crossover compounds with photo- 
magnetization, proposing a new model in which the spin-crossover phenomena 
and magnetic ordering are incorporated in a unified way. By using this model, 
the novel characteristics observed in Prussian blue analogs are qualitatively 
well reproduced. We investigate the time evolution of the magnetization and 
high spin fraction taking into account multi-time scales in a master equa- 
tion approach for the magnetic relaxation, the lattice (electronic) relaxation, 
and the photoexcitation process. In particular, processes with different tem- 
perature cycles starting from the photoinduced saturated magnetic state are 
studied including the effect of photoirradiation. In the temperature cycle in 
the low temperature region where the high spin state has a strong metasta- 
bility, the magnetization almost exactly follows the quasi-static process. On 
the other hand, when the temperature is raised above the region, the dy- 



namics of the spin state and the magnetization couple and exhibit various 
types of dynamical cooperative phenomena under time-dependent control of 
temperature and photoirradiation. 
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I. INTRODUCTION 



Spin-crossover (SC) transitions between the low spin (LS) state and the high spin (HS) 
state are induced by changes in temperature, pressure, photoirradiation, etc. [1] because the 
LS or HS state is realized due to subtle balance of the interaction between the electrons 
and the ligand field. Phenomena originating from these transitions have attracted much 
attention because they would potentially produce various functional devices [2]. Among 
them, photoinduced phase transition is a current topic in the field of molecular metal complex 
solids. There, much effort has been dedicated to reversible controls of magnetic, electronic, 
structural properties, etc. [1-9]. Molecular SC solids and Prussian Blue analogs (PBA) are 
typical materials that exhibit such properties. In order to understand the phenomena, it 
is necessary to study the dynamical processes of the transitions of photoexcitation and the 
relaxation of electronic states, where various cooperative effects play important roles. Thus, 
photoinduced non-equilibrium dynamics is now an area of intense study [10-17]. 

It has been clarified that SC transitions are not assembly of independent transitions 
between two electronic states in individual molecules but cooperative phenomena in cor- 
related molecules between a nonmagnetic (LS, diamagnetic) phase and a magnetic (HS, 
paramagnetic) phase [1]. Wajnflasz and Pick (WP) [18,19] theoretically gave the ground 
of the mechanism of the cooperative transitions by using an Ising model with degeneracies. 
The WP-model and extended WP-models have explained successfully static and dynamical 
properties of various cases of SC transitions [20-25] . 

The discovery of LIESST (lightinduced excited spin state trapping) [1,3] has led to a 
number of studies concerning relaxation during and after irradiation. The self-accelerated 
relaxation of HS fraction after irradiation and lightinduced thermal hysteresis, which are 
typical features in experiments, were reproduced by an Ising-like model based on the WP- 
model with Arrhenius dynamics [21]. To reproduce the relaxation dynamics with a slowing 
down at long time, a treatment beyond mean-filed approaches is necessary [24,26,27], where 
short range correlations (and sometimes long range correlations) between molecules and 
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spatially non-uniform fluctuation are important. 

In a series of related experiments, a photoinduced magnetic phase transition was ob- 
served in Co-Fe Prussian blue analogues (PBA) [4,5]. At low temperatures, PBA shows the 
transition: Fe n (S=0) + Co m (S=0) -> Fe n (S=l/2) + Co m (S=3/2) by photoirradiation. In 
the equilibrium, the system is in the LS (the left hand side) state which is a nonmagnetic 
state. By the photoirradiation, the system has a spontaneous magnetization where the spin 
state is in the HS (the right hand side) state. Furthermore a reversible switching between 
the LS and HS states by irradiations with different frequencies was observed in this com- 
pound. There, the magnetic interaction between the high spins should play an important 
role in addition to the lattice interaction which drives the LS-HS transition. Then, to study 
such multi-functional phenomena, it has become indispensable to establish a unified formu- 
lation which allows us to treat both SC transition and magnetic order-disorder transition 
simultaneously. 

In order to investigate the reversible switching mechanism of magnetization, the Blume- 
Capel (BC) model [28,29], which has a long-lived metastable ferromagnetic state in a param- 
eter region, was studied [30,31]. There, the dynamical process of the two-way switching was 
explained by focusing on the magnetic interaction and the multi-time scales of microscopic 
processes, e.g., the LS-HS transition, spin flip, and photoexcitation. 

For cooperated phenomena of the spin- crossover transition and magnetic ordering, we 
have proposed a unified model in which both the inter-molecular lattice interaction and the 
inter-molecular magnetic interaction are taken into account [32] . This generalized model in- 
cludes both the WP and BC models as special cases and enables us to treat the spin- crossover 
phenomena and magnetic ordering in a unified way. We have studied static properties of the 
model, and obtained its phase diagram as a function of the magnetic and lattice coupling 
parameters [32]. 

In this paper, we study the dynamical properties of the generalized model. In particular, 
we investigate the dynamical properties of the magnetization and the HS fraction under tem- 
perature cycles after photoirradiation, and we also study the process under photoirradiation, 
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motivated by the available experiments on PBA [11-13]. 

In order to study the dynamical aspects, we consider a master equation with multi-time 
scales where the relative relation of the relaxation times of the processes, e.g., transition 
between the LS and HS states, spin flip, and photoexcitation, plays an important role. We 
adopt here Arrhenius dynamics to realize the slow dynamics of the LS^HS transition at low 
temperatures, while we use standard Glauber dynamics for fast relaxation of the spin flip 
process. The process of photoexcitation is taken into account by a local excitation process 
from the LS state (0) to the HS state with up or down magnetization (±1) with the same 
probability. The strength of the irradiation is controlled in our formulation. 

In the present analysis, the following properties are clarified. Because the interaction 
between spins is considered, the typical shapes of the increase of not only the HS fraction 
but also magnetization by photoirradiation at low temperatures are naturally observed as 
non-linear functions of irradiation time [11,33]. Experimentally, it has been observed that 
the HS fraction and magnetization change drastically depending on the temperature cycling 
processes [11-13]. Here, we theoretically study the mechanism of the dynamics under various 
temperature cycles. In a cycle, we first prepare saturated magnetization under photoirradia- 
tion, and after turning off the irradiation the system is warmed up to a certain temperature, 
and it is cooled down to the initial temperature. It is found that in the low temperature re- 
gion where the excited high spin state has a strong metastability, the magnetization almost 
exactly follows the quasi-static process. On the other hand, if the temperature is raised 
above the region, the dynamics of the spin state and the magnetization couple and exhibit 
various types of dynamical cooperative phenomena under time-dependent control of tem- 
perature and photoirradiation. Our model explains qualitatively the novel behavior found 
in the experiments [11-13], and also predicts new phenomena which give useful suggestions 
for related experiments of photoinduced hysteresis phenomena [7,8]. 

The layout of this paper is as follows. In Sec. II, the model and method used in this 
study are explained. In Sec. Ill, dynamics in various temperature cycles are presented. 
The dynamical process under photoirradiation is given in Sec. IV. Section V is devoted to 
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summary and discussion. 



II. MODEL AND METHOD 

We have to introduce a model with both the magnetic and elastic interactions for in- 
vestigation of both SC transition and magnetic ordering. The degeneracies of the HS and 
LS states play an important role in SC transitions [18,19]. The effect of the degeneracies is 
expressed by a temperature dependent field [21]. In the present work, we adopt the following 
model [32,34]: 

n = -J £ sft - k E(^ 2 - his 2 , -h + (D- 1 ^) £0S? -h-h^Si, S l = 0, 1, -1, 

(id) ^ ^ Pi i 

(1) 

where denotes nonmagnetic state (LS state) and 1 and —1 denote up spin and down 
spin states (HS state), respectively. Here, J is the magnetic interaction, D{> 0) is the 
ligand-field energy, g is the effective degeneracy ratio between the HS and LS states [21], 
h is the magnetic field, and (5 = Although real PBA has a sub- lattice structure and 
a ferrimagnetic state, we simplified it to a simple cubic lattice and adopted ferromagnetic 
coupling (J > 0) with the nearest-neighbor interaction between molecules (sites). The 
inter-molecular interaction originating from the elastic interaction is denoted by K{> 0). 
Here, we adopt an attractive interaction for the same state (LS-LS or HS-HS). We set the 
parameters as J = 5K, K = 300K, D = 1200K which give realistic values of transition 
temperatures for PBA compounds [11-13]. The effective degeneracy ratio g is estimated 
to be 100~1000, which is attributed to the degree of freedom of spin, orbital, and lattice 
and molecular vibrations, and is also temperature dependent [21]. Here for simplicity, we 
consider g as a constant and set g = 200. 

To describe the dynamics of the model, we adopt the stochastic dynamics described by 
the following master equation: 

j t P(S 1 ---S k ---S N ,t) = L(t)P(S 1 ---S k ---S N ,t), (2) 
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where P is the probability of the state {Si • ■ ■ S k • • • S N } at time t, and L(t) denotes the time 
evolution operator which consists of the transition probability per unit time Wk(Sk — > S' k ) 
at the fc-th site, where S k (S' k ) takes a value , 1, or — 1. 

In order to explain the dynamics of PBA, we have to introduce multi time scales for 
different relaxations [31]. The time scale of spin inversion (S = ±1) is considered to be 
much faster than that between LS and HS transition (S = and \S\ = 1) because the latter 
is accompanied by a structural change. To describe the relaxation process between the HS 
and LS states adequately, we adopt Arrhenius dynamics [21]. On the other hand, we adopt 
standard Glauber type evolution between S = 1 and S — —1 for fast relaxation of the spin 
flip process. With these choices the difference between the dynamics of LS^HS transition 
and spin flip can be effectively taken into account. 

In Arrhenius dynamics, the transition probability is given in the form: 

Wf*(S k - S' k ) = \ exp(-(3(E - E(S k ))), (3) 

where Sk is (±1) and S' k is ±1 (0), and E(Sk) is a local energy of eq. (1). The quantity 
E — E(S k ) represents the effective energy barrier between initial (Sk) and final (S' k ) states 
and Eq is set as E = 1000K [21,24]. We denote the time evolution operator L(t) of this 
process by L x&tt . 

The transition probability of Glauber dynamics is given by 

Wr(S k -> S' k ) = \(l - S fc tanh/3A£), (4) 

where S k = 1 (—1) and S' k = —1(1) and AE = J J2 Si- We denote this process by L spin . 

From the experimental data [11], we assume that the photoirradiation causes a local 
excitation of LS — > HS, i.e., S k = — > S' k = ±1 with the same probability: 

W» hoto (0 ±1) = \. (5) 

We denote this process by L p hoto- 

The total time evolution including all processes is given by 
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L(t) — Piatt L\a,tt(t) + Pspin-^spin(^) + Pphoto photo 

(*), (6) 

where the coefficients are related to the relative strength of the processes. Here, p p hoto is 
related to the intensity of irradiation. We control the intensity of irradiation by j9 p hoto under 
the condition: p latt + J9 spin + p p hoto = 1- In this study, we take Piatt/Pspin = 1- Other choices 
such as piatt/Pspin = 0.5,2, etc. do not change the qualitative properties of the relaxation 
dynamics, because the difference of the time scale is taken into account in the choices of 
L latt (t) and L spin (t). 

The order parameters of the present model are the magnetization (0 < m < 1) 

m = M> (7) 
and the high spin fraction (0 < q < 1) 

9 = %^, (8) 

where iV is the system size (set as iV = 20 3 ) and (• • ■) represents an average over the 
stochastic process. 



III. DYNAMICS IN THE TEMPERATURE CYCLES. 



In this section, we study three kinds of temperature cycles (heating and cooling) after 
switching off irradiation to simulate the experiments [11,12]. First, the system is irradiated 
(Pphoto = 0.01) at a very low temperature (T = 5K) under a magnetic field (h = 0.1) until 
the state becomes the saturated magnetic HS state (q — 1 and m — 1). 
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FIGURES 

FIG. 1. Time dependence of m (solid line) and q (dashed line) in the initial photoexcitation at 
T =5K. 

In Fig.l we show the time evolution of m (solid line) and that of q (dashed line) at 
T under irradiation with h = 0.1. The change of m follows that of q because at this low 
temperature the magnetic moments of excited HS states rapidly align in the direction of the 
magnetic field due to the magnetic interaction. Here, we find characteristic time evolution 
of m and q. That is, at the beginning of the irradiation, m and q increase slowly with 
concave shapes. Then, they are accelerated and finally saturate slowly with convex shapes. 
This time dependence is consistent with the experiments [11,33]. This result indicates the 
following. At the beginning of the photoirradiation, the irradiation makes local excited HS 
states but the cooperative effect tends to put them back to LS states because almost all 
neighboring sites are in the LS state, and thus the growth of q and m is suppressed. It 
should be noted that this initial suppression is due to the cooperative interaction between 
LS sites. In the present scheme, this effect is automatically taken into account. When the 
number of HS sites is increased, the cooperative effect between HS sites becomes significant, 
resulting in the acceleration of q and m. 

After switching off the irradiation and field (setting p p hoto = and h = 0), we change the 
temperature. The temperature is gradually increased up to T max , and then decreased to the 
initial temperature T (i.e., T — > T max — > T ). Here we change the temperature in IK steps. 
At each temperature, first 1000 Monte Carlo steps (MCS) are discarded for a transient time, 
and then 1000MCS is used to measure the quantities. Thus, the temperature dependence 
of the quantities denotes their time dependence. The following three cases are studied. (A) 
T max = 60K, (B) T max = 250K and (C) T max = 93K. We will see that, depending on the 
value of T max , the system shows qualitatively different dynamical behavior of q and m. 

FIG. 2. Time (temperature) dependence of m (solid line) and q (dashed line) in process A. The 
arrows with solid heads (bare heads) indicate the direction of change of m (q) . 



9 



The time (and temperature) dependence of m (solid line) and that of q (dashed line) 
in process A are depicted in Fig. 2. The magnetization m decreases with heating and 
disappears around 23K. The magnetization in the cooling process almost exactly follows 
that in the warming process. This indicates that the dynamics of m follows the quasi- 
static magnetization process. It is noted that g(~ 1) is constant during the whole process. 
Although the free energy in the static mean-field theory for the equilibrium state as a 
function of m and q does not have metastability for q = 1, the strong metastability of 
q — 1 appears as a result of Arrhenius kinetics with a relatively large E as mentioned in 
the previous section. Indeed, it is found that the metastability appears in the dynamical 
mean-field theory, where the time evolution of the magnetization is given by the Van Hove 
equation: 



where a double well potential of m is realized. Thus, under the metastability of q(—l), m 
behaves quasi-statically and shows a typical temperature dependence of the second-order 
magnetic transition in both warming and cooling processes. Assuming q — 1, the transition 
temperature is given by the Curie temperature of 3D Ising model on the simple cubic lattice 
[35]: Tq = 4.5J = 22. 5K. This quasi-static behavior is consistent with experiments [11,12]. 

FIG. 3. Time (temperature) dependence of m (solid line) and q (dashed line) in process B. 

Next, we study a case with a high value of T max (T max = 250K, case B). Figure 3 shows 
the time (and temperature) dependence of m and q. Now we find self-accelerated (sigmoidal) 
relaxation of q in the region 90-110K, above which q disappears until the SC transition occurs 
at a higher temperature. We find hysteresis behavior in the region 160-220K which is due 
to the SC transition. The dependence of q and that of m on the temperature cycle have 



m = 



df(m,q) 



(9) 



dm 



As far as q ~ 1 at low temperatures, the free energy is given by 




(10) 
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been found in the experiments of PBA [11-13]. In the process, the magnetic interaction has 
negligible contribution to the thermal transition of q because of J <C K. Thus the transition 
point of the spin crossover where q — 1/2 is given by 



exhibits the first-order phase transition at T cross [18-24]. Thus, the center of the hysteresis 
depends on D and g and the width of the hysteresis depends on the ratio K/D for fixed g 
and the changing rate of the temperature, which indicates varieties of hysteresis loops in the 
experiments [13]. After arriving at T max , the system is cooled down. In the cooling process 
below the spontaneous thermal hysteresis region, m and q are zero, where the system is in 
the equilibrium state. 

FIG. 4. Relaxation of q when the temperature raising is stopped at several temperatures (80, 
85, 90, 95K) in process B. 

If we stop warming at some temperature on the way to T max from T in process B and 
maintain the temperature, we observe characteristic sigmoidal relaxations of q with a tail as 
depicted in Fig. 4. These relaxation curves with temperature dependence are also observed 
commonly in the related experiments [7,9,21]. 

FIG. 5. Time (temperature) dependence of m (solid line) and q (dashed line) in process C. 

Let us consider the third case C (T max = 93K). Just before arriving at this T max in the 
warming process, q begins to relax. In Fig. 5, the time (and temperature) dependence of m 
(solid line) and that of q (dashed line) are plotted. The relaxation of q is observable until T ~ 
60K in the cooling process. The quantitative relaxation profile depends on the changing rate 
of the temperature, but qualitative features are common. When the temperature is further 
decreased, q freezes with a non-saturated value of q(< 1). Because q < 1, the substrate 
for the magnetic order is in a sense diluted, and the magnetic transition temperature in 
the cooling process is lower than the Curie temperature Tc for q = 1. In this process, the 




(11) 



Because 



Across is less than the Ising phase transition point: 4.5 > T cross / 



j, the system 
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fluctuations of m and q are large. Thus, we take average over ten independent simulations, 
and plot the average with error bars in Fig. 5. This lowering of the magnetic transition 
temperatures was observed also in the experiments [11,12]. 

IV. DYNAMICS IN THE TEMPERATURE CYCLE UNDER 
PHOTOIRRADIATION 

FIG. 6. Effect of the irradiation (p p hoto = 0.0001) on m and q in process B. (a) q under 
continuous irradiation (dotted line) is compared with q with no irradiation (dashed line). The 
arrows indicate the direction of change of q in the former case, (b) m (solid line) and q (dotted 
line) in the low temperature region under continuous irradiation. 

In this section, we investigate the time (and temperature) dependence of the order pa- 
rameters under photoirradiation, which would give useful suggestions to related experiments. 
We study the same temperature cycle as B under photoirradiation (p p hoto = 0.0001). In Fig. 6 
(a), the time (and temperature) dependence of q under the irradiation (dotted line) is shown 
in comparison with the data in the case of no irradiation (same dashed line as in Fig. 3). Be- 
cause photoexcitation enhances the HS state and suppresses the relaxation of q, the thermal 
hysteresis of q in the region 160-220 K is shifted to the low temperature side. In addition, 
lightinduced hysteresis is observed at lower temperatures (10-120 K). The sigmoidal relax- 
ation of q, which is the right branch of this hysteresis, starts at a higher temperature than 
in the case with no irradiation. This tendency is consistent with the dynamical mean-field 
analysis of the Ising-like model [21]. 

In Fig. 6 (b), the time (and temperature) dependence of m (solid line) and that of q 
(dotted line) are shown, focusing on the low temperature region. In the warming process, 
m disappears at a transition temperature which is the same as To Indeed we find eq. (10) 
has no p p h to-dependence when q — 1 under irradiation (p p h to 7^ 0). It is worth noting that 
in the cooling process, m appears at a temperature lower than Tc(— 4. 5 J = 22. 5K). The 
irradiation causes growth of the HS fraction q around 30K and after a delay the magnetic 
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order appears at T' c ~ 18K. Before the full saturation of q, the HS sites have a diluted 
distribution (q < 1). Thus, the critical temperature of the magnetization is lower than T c . 
At the critical point T c , the high spin fraction is q ~ 0.6 in Fig. 6 (b). 

For this value of q, the magnetic order should appear at 0.53T C [36] in the randomly 
diluted lattice (percolated lattice). The present T c (~ 0.8T C ) is higher than this value, 
which should be attributed to correlated high spin states. The nature of this correlation of 
the HS sites will be studied in detail elsewhere. 

The present observation suggests that if the temperature at which the left branch of the 
photoinduced hysteresis of q appears is near Tc, we find the magnetic hysteresis driven by 
the photoinduced hysteresis of q in the temperature cycle. If q appears at a much higher 
temperature, the temperature dependence of m is the same as that in the warming process. 
We expect that this magnetic hysteresis will be found in experiments by controlling the 
strength of irradiation. 

V. SUMMARY AND DISCUSSION 

We have presented a three-state model which allows us to describe both SC and magnetic 
properties of PBA with photoirradiation. By using this model, the time (and temperature) 
dependence of magnetization (m) and that of the high spin fraction (g) were observed 
under several temperature cycling processes. In the case of the temperature cycle in the 
low temperature region (T max = 60K), the second order magnetic transition takes place 
at the Curie temperature (T c ~ 22. 5K) of the corresponding 3D Ising model in both the 
warming and cooling processes, because the strong metastable high spin state is realized 
(q ~ 1) in this temperature region. When T max is high (250K), q disappears showing the 
self-accelerated (sigmoidal) relaxation, and it appears again at a high temperature and 
the thermal hysteresis of spin crossover is observed. In the cooling process, the system 
is in equilibrium (diamagnetic state) below the thermal hysteresis region. When T max is 
moderate (93K), partial relaxation of q is observed in the cooling process, and a magnetic 
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phase transition occurs at a temperature lower than Tq- There, the saturated m is less than 
the fully saturated m{= 1) because the saturated q is less than 1. Our results reproduce 
the experimental results [11,12] qualitatively, and thus we believe that the dynamics that 
we adopted here represents an important aspect of the real dynamics of PBA. 

We have also investigated the dynamical properties in the temperature cycle under ir- 
radiation. A shift of the sigmoidal relaxation of q to a higher temperature and a shift of 
thermal hysteresis to a lower temperature are observed, because the photoexcitation makes 
high spin state more stable. We found that the magnetic hysteresis caused by the photoin- 
duced hysteresis of q occurs under a special condition, which is expected to be observed in 
related experiments in the future. 

Very slow relaxation phenomena in photoinduced magnetization, such as the different 
values of magnetization in ZFC and FC procedures like the spin-glass systems, have been 
reported in some PBA very recently [11,12]. With the present model, those phenomena will 
be analyzed, which will be reported elsewhere. 
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